IP-ASTP-3-97 



April, 1997 



Double parton scattering of hadron-hadron interaction 
and its gluonic contribution 



Hung Hsiang Liu [] 

Institute of Physics, Academia Sinica 
Nankang, Taipei, Taiwan 115, R.O.C. 



Abstract 



We propose a formalism for calculating the event cross-section of parton interactions. In this 
formalism, we use the light-front bound-state wave function to expand scattering initial state. This 
leads to an expression of the cross-section in terms of square of wave functions which can be used to 
define the parton distribution functions. The probability interpretation of the partonic interactions 
among hadron-hadron scattering is naturally achieved. We apply this formalism to calculate the 4- 
particle final state cross-section for single parton, double parton and gluon splitting interactions. We 
compare the behaviors of these cross sections and propose a method for experiments to differentiate 
these interactions. 
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I. Introduction 



Phenomenological calculation of hadron-hadron cross section uses the probability of 
finding one parton in each scattering hadron multiplied by the parton-level cross section da: 

da(A + B -> n jets + anything) = ^ / dx a dx b F a/A (x a )F b/B (x b ) 

a,b 

x da (a + b — ► n jets). (1) 

In the calculation, one uses universal parton distribution functions for F and uses pertur- 
bative QCD to calculate the n-jet single parton cross section &2^n- Since 1982 various 
experiments studied two-, three- and four-jet events from hadron-hadron scattering Jl|. The 
two- and three-jet cross sections agree well with theoretical calculations ||. However, the 
experiments on four-jet events show some inconsistence to the calculations ||. The UA2 M 
experiment showed a good agreement to the calculations based on single parton interaction. 
On the other hand the pp collision experiments by AFS Collaboration at CERN ISR || 
and the CDF experiments |J showed a small violation to the calculations. Consequently 
CDF and AFS collaborators suggested that a small double parton contribution to the QCD 
(2 — ► Ajets) cross section may provide a better fit to the experimental data. 



The double parton interaction happens when two partons in one hadron collide with 
two partons in another hadron [Fig. 1]. Since this process involves more than one parton 
in a hadron at one time, it will give more informations on hadronic structure than the 
single parton scattering can give. Experimentally one can differentiate the double parton 
events from single parton events by the final jet configurations. First, since the transverse 
momenta of partons in a hadron are small, the transverse momenta of the two jets produced 
in each of the parton-parton scattering should be well balanced. Therefore in a double parton 
scattering, one can find two pairs of jets with transverse momenta, say (pi, P2) and (P3, P4), 
such that |pi + P2I, |p3 + P4I <<1 GeV. In the contrast, in an ordinary four-jet event one 
will have Pi + P2 + P3 + P4 ~0, but for any pair of transverse momenta (Pi,Pj), one will 
normally find |pi + Pj| >>1 GeV. Secondly, for double parton scattering, if we assume that 
there should not be strong correlation between the two partons in the individual colliding 
hadron, the distribution of the double parton cross section in the azimuthal angle between 
the jet pairs should be flat. The CDF collaboration defines two variables which are sensitive 
to these signals. They define variable S to exploit the tendency of pairwise balance, 



S(i + j,k + 1) = min 



k 'V^' p )2+( / l IP , 2 +Pd ^ (2) 
Y viPii 2 +iPji 2 vi p *i +Ip*i 
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where S is minimized over the three possible jet pairings (12,34), (13,24) and (14,23) and 
they define the variable AS to be the azimuthal angle between the two jet pairs which 
minimize S. The double parton events should favor smaller values of S and the single 
parton events should favor larger ones. Also, the double parton events should be uniformly 
distributed over AS in the interval to n. 

In this paper we discuss the theory of the 4-particle final state events from hadron- 
hadron scattering. Most of these events are in the form of 4-jet production. There is 
a small fraction of these events which is in the form of di-lepton pair production. It is 
believed that the hadron fragmentation of the jet final state, which is a soft process, does 
not change the kinematics of the final particles. Therefore the di-lepton pair production 
should have a similar kinematics to the jet production. For a theoretical understanding of 
double parton scattering, we will use the di-lepton process as examples to the formalism. 
The same formalism can be applied to calculate the jet cross section. We are interested in 
the experiments with observables Qi, Q 2 of the TeV order, where Qi and Q 2 are the invariant 
momenta of the final particle pairs. The largest source of these events is the lowest order 
single parton scattering. Double parton scattering may contribute these events. There is 
also a possible type of four-particle events that a gluon splits before collision. In this process, 
one incoming gluon splits into two partons, while the other incoming gluon also splits, then 
the two pairs of partons collide. One can interpret this process as an alternative form of 
double parton scattering. Imagining that the parton pairs are produced in time At before 
the collision, where At is long compared to the parton collision time 1/Q but short compared 
to the typical time for partons within a hadron to interact, then the transverse momenta of 
the partons produced in the gluon splitting process will be in the rang of 0<< k << Q i: Q 2 . 
We will examine the effect of these large transverse momentum partons on double parton 
scattering. The paper is organized in the following fashion. In the next section, we introduce 
the light-front wave functions which are useful for the formalism. In section III, we derive 
the cross section for the single parton scattering. In section IV, we use the formalism 
to calculate the double parton cross section. In section V, we discuss the gluon splitting 
scattering. Finally, in the discussion section we will compare the four-particle final state 
cross sections of single parton scattering, double parton scattering and from gluon splitting 
scattering. We propose a method for experiments to differentiate these interactions. 
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II. Light-front wave functions 



In the light-front field theory, one uses a Hamiltonian approach to derive the equation 
for a bound-state wave function. For this purpose, one uses gauge field theory quantized 
on light-front planes. That is, one attempts to solve the eigenvalue equation P op | x I / ) = 
E\^}, where P Q ~ is the generator of translations from one light-front plane x + = const 
to the next. Consider a state |P + ,P) which denotes a hadron with light-front momentum 
(P+, (P 2 +M 2 ) /(2P+), P), normalized by (P+, P|P'+, P') = (27^) 3 2P+5(P+-P'+)5 2 (P-P , ). 
We can expand the state vector in the Fock space 

\P + ^) = E E (f[^^)b + (pt,PuSi)..d + (p+,p m ;s m )\0) 

m m 

x -p m , si, ... Sm )2P + (2vr) 3 5(P+ - $> t +)5 2 (P " E Pi)- ( 3 ) 

i=l i=l 

The advantage of this framework that one hopes will make the problem tractable is that 
partons (at least those with zero momentum) cannot appear spontaneously from the vacuum: 
The internal state of the hadron is given in terms of a set of wave functions ip n defined by 

(0\b(pt, Pi; s!)d(pt, p 2 ; s 2 )...d(p+, p m ; s m )\P + , P) 

m m 

= (2tt) 3 2P+5(P + - $>+)5 2 (P ~ E Pt)Mpt, Pi;P 2 + , P2-P+, p m ; «i, ..s m ), (4) 

i=l i=l 

where 6's and cPs are the destruction operators for quarks and antiquarks. The indices Sj in 
the wave functions are the internal symmetry variables which denote helicity A«, isospin 
and color q. From the normalization equation we can show 

i = E E (n7^^)(2-) 3 2P + 5(p + -E^ 2 (P-EpO 

m si...s m i=l Z A i=l i=l 

X |^m(Pl",Pl,--Pm,Pm;Sl,...S m )| 2 . (5) 



The normal light-front wave function is good for describing particles moving in the +z 
direction. In order to describe the bound-state particles moving in the — z direction, we have 
to do the followings: We expand the state vector for hadron moving in the — z direction in 
the Fock space by 

l p ~ p > = EE(n^# i ) &+ (pr,Pi;«i)--rf + (p;,Pn;«n)|o) 

n n 

x Vn(pi, ..Pn\ si, ...s n )2P-(2nf5(p- - J2PI )S 2 (P - E P*)- ( 6 ) 

1=1 i=l 
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Define the spinors fi(p, s) and v(p, s) as the spin components of quarks and anti-quarks 
moving in the —z direction: 



//(p, 1/2) 



/ p 1 — ip 2 \ 
V2p- 


\ m / 



/x(p,-l/2) 



/ m \ 


V — p 1 — ip 2 / 



v(p, 1/2) 



/ — m \ 


\ — p 1 — ip 2 / 



f(p, -1/2) 



\/2p- 



( p 1 — ip 2 \ 
V2p- 


V — m 



(7) 



The normalization of these vectors is: 

J2Va(p~,P,s)jI a (p-,PiS) = -J2v a (p~,p,s)v a (p~,p,s) = j4 + m. (8) 

s s 

It can be shown that these spinors are the eigen-states of a helicity operator referring to a 
Lorentz frame moving in the +z direction near the speed of light, i.e., v z = — tanh(cu), uo — > 
— oo |7]]. The helicity operator is 



/ P 1 - W 2 



Kip) -> h-ooip) 



V2p- 






V o 











o \ 








-p~ 







+ ip 2 


V— J 



(9) 



One can check that the spinors listed in Eq. (0) are eigen states of h_ OCJ ip): 

h-ao(p)v(p, = T^(p, ±~)- 



(10) 



We therefore clarify that the physical meaning of the spin index s is the helicity of the quark 
as measured in a Lorentz frame moving near the speed of light in the +z direction. Similarly, 
one can define the photon polarization vectors e^(A) as 



The physical meaning of the index A is the helicity of the photon as measured in a frame 
moving the +z direction. 
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III. The cross section for single parton scattering 



In this section we derive the cross section formula for single parton scattering. For 
simplicity, we use the di-lepton production from unpolarized hadron-hadron scattering as an 
example. We are interested in the experiments with observables Q\, Q2 of the TeV order. 
The transverse momenta of partons in a hadron are therefore negligible. The differential 
cross section for this process as shown in Fig. 2 is 

m n 

x (2vr) 4 (5 4 (P A + P i? -^a l -^6 J -g 1 -g 2 ), (12) 

i=2 i=2 

where Q\ and Q2 are the invariant 4-momenta of the lepton pairs I2) and U) and 
|M| is the invariant hadronic scattering amplitude which is calculated by using its relation 
to the scattering matrix 

out (anything, l u l 2 , l 3 , h\S\A, B)- m = (2n) 4 5 4 { Pi - Pf )M. (13) 

The state \A, B)- in denotes the initial state for hadron A moving in the +z direction and 
hadron B moving in the — z direction. We assume that the state \A,B) in can be factor- 
ized into \A) in \B) in . Then we expand the hadron state \A)i n by Eq. ([|) and expand the 
hadron state \B)i n by Eq. (J6|) After applying these expansions we can express the scattering 
amplitude square as 

m r rln + rl 2 a- n r rlh~rl 2 h- 

mi = eeh En/ ^Ifi-i 2 

m,ns a2 ..s an i=2 J l Z7i J ZU i s b2 ..s bn i=2 J ^ Z7t ) LV % 

^ 2~ ^Pam ("^al 5 ^ ] Sali ^a2i 3-2; ^a2---j •^•ami Sam) 
X al i=2 

m 

X — i>bni X bl, Ql + Q2 + a " S 6l! X bn, b 2, S b2 ...; b n , S bn ) 
X bl i=2 

x /d 2 Pa 5 2 ( Pa -f>). (14) 

J i=2 

The wave function ipami^bn) is a function of the hadron spin Sa(Sb), but we do not explicitly 
display this dependence. We inserted an identity 1 = / d 2 p a 5 2 (p a — X)£L 2 a «) m Eq. ( |14"D for 
future purpose. This insertion serves to define p a as the total transverse momentum of the 
non-interacting partons numbered from 2 to m. The parton momentum fractions are defined 
by x ai = a\ jP\ and x bi = /Pq. The momentum fractions of the interacting partons ai 
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and b\ are fixed by the overall delta function. They are 

M\ hi. 



pi 


a t 


Qi + Qi 




pi 


p + 


Pb 


- £?= 2 K 


QI+Q2 




Pb 


Pb 



x al = a ^:= z 1 = +0 

s s 
M 2 a 2 

x bl = ~ a Z- = ' Z " % = ^ l n ^' z +0(^,^). (15) 

In this equation |m| 2 p is the invariant partonic scattering amplitude square. Its integration 
over the final particle phase space denoted by |T| 2 p will be calculated in the appendix. In 
order to simplify the formula, we collectively express the sets of the non-interacting parton 
momenta and spins by: 



1=2 



i=2 



Va 


= {(4 


, a 2! •• 


■ d m , a m 


>IE^ 

i=2 


v b 


= {(b 2 , 


b 2 ; .. 


■K., h n) 


n 
i=2 


S a 


— { s a2, 


• ■■Sam 


m 

. | ^ t $ai 
i=l 


= S A } 


s b 


= {Sb2, 


■ ■■Sbn | 


n 

E^ = 


= S B }, 



i=l 

and we denote the phase space of the non-interacting final partons as 

r m r rln + rl 2 a m m 

/ W a ] = n / ^4^ 2p A^(Pi -H a t-Qt- Qt)s 2 (p a - £ * 

J i=2 J \ Z7V ) Za i i=2 i=2 



(2vr)32a+ 

J [dVb] = tlJ§^ 2P B( 2 ^ 6 ( p B-tK-Qi-Q2 



With these notations we can rewrite Eq. (|r2|) as 
da da 



d 4 Q 1 d 4 Q 2 d 2 Q T d 2 AQdQfdQidQtdQ2 

1 1 ^E^EE/^K]/^] 



(16) 



x 5 2 (Pa + Qi + Q 2 + £*>*). ( 17 ) 

i=2 



4s 2 (27i) 2 x 2 al xl2 
x E E \ T \ 2 [ d2 Pa /dVPfo. + ai ) 

m,ns al ,s bl 

x ^lm( x ai, &u s al ;p a ;V a , S a ) 

x ipl n (x al , Q T - ai, s bl ; -Q T + p a ; V b , S b ), (18) 

in which we define Qt = Qi + Q29 AQ = (Qi — Q 2 )/2. For the following discussion we 
inserted an identity 1 = / <i 2 ai5 2 (ai + p a ) in Eq. ([TJ) The wave function ip am (or ipbn) in 
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Eq. ( |T8D is the same wave function in Eq. (|T4| ) expressed by the new variables ai, p a , V a (or 
bi, pt,, V b ). We can use these momenta to define the following variables: 

Ap a = x a a 1 - x al p a 

A = p a + ai, (19) 

in which x a is the total momentum fraction of the non-interacting partons and A is the total 
transverse momentum of hadron A. By using coordinate r tt x of parton a\ and the center 
mass coordinate r a of the remaining partons, we can define the following variables which 
conjugate to momenta Ap a and A: 

Ar a = r a x — r a , 

R a = x al r al + x a v a , (20) 

here Ar a is the relative transverse distance between parton a\ and the center mass of the 
non-interacting partons and R a is the center mass coordinate of hadron A. Similarly, we can 
define the momentum variables Apf, and B and their conjugate coordinates Ar b and for 
hadron B. To express the cross section Eq. ( |T8| ) in terms of these coordinate variables, we 
first express the wave functions in terms of Ap a and A then Fourier transform these wave 
functions to the coordinate space by using, for example, 
V>(Ap ,A;F A ) = l/27r/rf 2 Ar a exp(zA Pa - Ar a )^(Ar a , A; V A ). 



The integration volume / d 2 p a J d 2 a.i in Eq. fll8|) can be replaced by / <i 2 p a 2i J d 2 A. It 
is easy to see from the overall delta function that the total transverse momentum Qt is of 
the order of the parton transverse momentum which is of a few MeV. Since we are interested 
in the experiments with observables Qf, Q~ of the TeV order, we can neglect the p 2i, A 
and Qt dependence in \T\ 2 and we can integrate Eq. (|18|) over these transverse momenta. 
The integration of Qt gives us a delta function which enables us to express the differential 
cross section in terms of the wave function squares: 

da 1 1 l_y^ ^ , , 2 

d 2 AQdQtdQ{dQUQ2 As 2 {2ix) 2 x 2 al x 2 bl ^ n! 2 ^ 1 1 



'. S a \S bl 



oEE \dV a ] / d 2 kr a - a $ 

am\X a li s a \, Ar aa , V a , S a ) I 

OA So 

^EE/ Wb] J d 2 Ar bb \^ bn (xbi,s bl ; Ar b - b ;V b , S b )\ 2 . (21) 



2 

Sg St 



We can define the following function f a (and similarly f b ): 



m S A Sa 
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X 1*0 am\%al} S a \ } Av aa ,V a , S a )\ . 

(22) 



From the normalization equation Eq. ([5]) , one sees that 

E/dW.(Wox) = l. 



(23) 



The function / a can then be interpreted as the distribution function for parton a\ with 
momentum fraction x a \ and internal symmetric variable s a \- 

If we assume that for non-polarized hadron the distribution functions of different helic- 
ities are equal, i.e., f(x, 1) = f(x, —1) = 1/2 f(x), we get the cross section for single parton 
scattering as 

da da 
d 2 AQdQfdQ^dQtdQ2 ~ \d 2 AQdQ 2 1 dy 1 \dQ 2 2 dy 2 

-f a (x al )f b (x bl )\T\ 2 , (24) 



(2vr) 4 1 1 ( ; , w , , ,,. r , 2 



4S 2 X a i X b i ' 

where y's are the rapidities defined by yi = 1/2 \a.{Qx/Qj) and \T\ 2 is calculated by using 
the (2 — > 4 particle) Feynman diagrams. 



IV. The cross section for double parton scattering 



In this section we derive the cross section formula for double parton scattering. For 
simplicity, we use the di-lepton production processes in unpolarized hadron-hadron scattering 
as an example. The process is shown as in Fig. 1. By expanding the initial states as in the 
previous section, we can express the scattering amplitude square for double parton interaction 

as 

m r dn + d 2 n- 11 r dh + d 2 h- 

mi = £ £ e n/g0 

m,ns a3 ..s an i=3 J J ZU J s b3 ..s bn i=3 J > LV i 



X (2nUs 2 x 2 U I ^ I 



X ^Pam^al: Sal) ^a2i a l ^ ] 3ij $0,2) ^a3) a 3; ^a3-- m ) ^ami ^mi ^am) 

i=3 

m 

X Ipbnfcbl, Ql - a l, ^62, Q2 + ai + Xl a *' S b2', ^63, b 3 , s 63 ...; x 6n , b n , s bn ) 



8=3 
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X :'a/rn ; a i> ^cslj X a2l <^\i / t 3tj ^a2' *^a3> a 3, ^a3 --- ) -^am; ^m; "^amj 

i=3 

m 

x ^bn( X bl, Qi - a 'l) 4l! ^62, Q2 + a 'i + E x 63, b 3 , s 63 ...; x„, b n , s 6n ). (25) 

i=3 

In this equation, |m|^ p is the partonic scattering matrix square and the integration over 
its final particle phase spaces, \T\^ , will be calculated in Appendix. To calculate the jet 
productions from double parton scattering, one can apply Eq. and use the corresponding 
|T| 2 for the jet production. The non-interacting partons are now numbered from i = 3 to 
m or n. The momentum fractions of the interacting partons a%, a 2 , bi, b 2 are again fixed by 
the overall delta function. They are 

x ai = Qf/P+, x M = Qr/P B , 2 = 1,2. (26) 



In order to simplify the expression for the cross section formula, we define analogically 
to the previous section the collective variables V a , V&, S a and Sb for the non-interacting 
partons. The only difference is that there are now m — 2 or n — 2 non-interacting partons. 
Therefore we can rewrite Eq. (|25| ) as 

da 1 1 

#Qi#Q 2 ~ 16s 4 (2vr) 6 x2 1 x 



2 „2 
62 



x ^ee^ee/to /wi e e m 2 

Sa S a Sb S b ' s a is a 2s' al s' a2 s b is b 2s' hl s' b2 

X J d 2 p a J (f&i J d 2 cL 2 5 2 (p a + & 1 + &2)4>a m (Qf,&l,Sal,Qt,&2,Sa2,Pa]V a ,S a ) 

x ipbn(QT, Qi - a i, s bl ; Q 2 , Q 2 - a 2 , s b2 , -p a - Qi - Q2; H, S&) 

x / d 2 p' a J d 2 a[ J d 2 a' 2 5 2 (p' a + &\ + a' 2 )^* am (Qf, a[, s' al ; Qt, a' 2 , s' a2 , p a , V a , S a ) 

x 4>L(Qii Qi - a 'i, 4i; Q2 7 Q2 - a' 2 , s', 62 -p a - Qi - Q 2 ; H, S b ) 

x 5 2 (Pa-p'J. (27) 



We inserted the following identities in Eq. (|27|) for a more symmetrical expression: 

1 = | d 2 p' a 5 2 (p' a - p a ) J d 2 a 2 5 2 (at + a 2 + p ) J d 2 a 2 S 2 (a' 1 + a' 2 + p'J. 
Define the following transverse momentum variables by using parton momenta a 1; a 2 and 



Pa- 



a 2 i — £ a ia 2 — £ a 2 a i> 

Pa21 = Xa{&l + a 2 ) - (X a i + X a2 )p a , 

A = p a + ai + a 2 , (28) 
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here we define £ a i = x a \j (x a i + x a2 ) and £ a 2 = $0,2/ (x a i + £02) ■ The variable x a is the sum of 
the momentum fractions of the partons which are not involving in the collision. Therefore 
we have £ a i + £ a 2 = 1, and x a \ + x a2 + x a = 1. The momentum A is the total transverse 
momentum of hadron A. By using the coordinate r tt i of parton a\, the coordinate r a2 of 
parton a 2 and the center mass coordinate r a of the remaining partons, we can define the 
following variables which conjugate to momenta a 21 ,p a2 i and A: 

?21a = I"o2 — I"al, 
'■aa "a 1-a, 

Ra = {x a i + x a2 )r a + x a r a , (29) 

here f a = £ a i r ai + £ a 2 r a2 is the center mass coordinate of the two partons in hadron A 
which are involving in the collision, r 21a is the relative transverse distance between these 
two partons and r Sa is the relative transverse distance between the center mass of these two 
partons and the center mass of the remaining partons. Note that these relative transverse 
vectors are Lorentz invariant quantities so that their values are the same as their values in 
a rest hadron. Finally, R a is the center mass coordinate of hadron A. The momentum and 
coordinate variables can then be replaced by the new variables by the following equations: 



a l — X a iA + £alP21a — a 21 , 




a 2 = x a2 A + £ a2 p 2 i a + a 2i , 




Pa = X a A - p 2 ia, 




r al — R-a + %a r aa ~ £a2 r 21aj 




Ta2 = Ra + %a r aa + £al r 21a) 




Ta Ra \%al Xa2)^aa- 


(30) 



The Jacobian of these replacements is 1. Similarly we can define the momentum variables 
b 2 i, P216, and B and their conjugate coordinate variablesr 2 if>, r& 6 , and R& for hadron B. To 
express Eq. ( P7D in terms of coordinate variables we Fourier transform the wave functions 
from the transverse momentum space to the coordinate space by using, for example 
-0(^21, P2i a ; ^4) = 1/ (2vr) 2 / c/r 2 i a (iraaexp(2a 2 i • r 21a + ip 2 i a • r go )V>(r 2 i a , r aa ; A). 
The integration volume / <i 2 p a / cP&i J d 2 a 2 in Eq. (p7f ) can be replaced by / d 2 p a2 i / d 2 a 2 i / d 2 A. 
Since we are interested in the experiments with observables Qf, Q~ of a few TeV, we can 
therefore neglect the p a 2i, a 21 , A dependence in |T| 2 . Integrating Eq. (|27D over these 
transverse momenta we get the differential cross section as: 

da da 
d A Qxd A Q 2 ~ d 2 Q T d 2 AQdQfdQidQtdQ 2 
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1 (Se4('0 7r \8 -.2 rj.2 „2 ™2 X! E l-^l 

1 1 f rnr 

X 



^ee^ee/wi/wi 

OA Oa OB Oft 

J d 2 v- aa J d 2 r 21a J d 2 r bb J d 2 r'- aa J d 2 r' 2la 



X 1pam(Xal, S a l] X a2 , S a2 , r aa ; r 2 la] V a , S a )lp bn (x b l, S bl ] X b2 , S b2 ] Y bb , Y 2la ] V b , S b ) 
X ^Pam( X ali s a l'i X a2, S a2 , T aa ; r 21a ] V a , S a )lp bn (x b l, S bl ] X b2 , S b2 , T bb ; r 21a ] V b , Sb) 

x exp{iQ T [r aa - r'- aa + -(£ lo - 6 )(r' 21a - r 21a )])exp[iAQ(r' 21a - r 21o )]. (31) 



It is again easy to see from the overall delta function that for double parton scattering 
process, both the momenta Qt and AQ are of the order of parton's transverse momenta 
which are of a few MeV. We can therefore neglect the Qt and AQ dependence in |T| 2 and 
we can integrate Eq. (|31f) over Q r and AQ. These integrations give two delta functions 
<5 2 (r' 21a — r 2 i a ) and 5 2 (r aa — r' aa ). We see from the Appendix that the spin flip process for 
|T| 2 p is negligible therefore we can set s a i = s' ai , su = s' M . These results enable us to express 
the differential cross section in terms of the square of wave functions: 

da = 1 1 1 J- \T\ 2 [d 2 r 21a 

dQ l dQ l dQ 2 dQ 2 16s 4 (27r) 4 x al x a2 x bl x b2 SalSa2SblSb2 J 

X 2 EE y ^ 2r aa|^ am\%al, S a l, X a2 , S a 2, ^aaj ^21a, 

V a ,S a ) 



2 



Sa S a 

1 

X 



7; EE [l dv b] I d 2 r bb \^ bn (x bl ,s bl ]x b2 ,s b2 ]r bb ]r 21a ]V b ,S b )\ 2 . (32) 

OB Of, 

We can define the following function 

fa{.X a i,X a2 ,Y 2 i a , S a i, S a2 ) = T 7o~^G 9 EE [ [^Va] f dr aa 

e ±X a \X a2 \£K ) Sa S a 

x \^am{^a\, Sal', x a2, S a 2) r aa', r 21a', Va, S a )\ 2 . (33) 

By the normalization equation Eq. (|5|) we can show that f a the distribution function for 
two partons with momentum fractions x a i and x a2 in a relative transverse distance r 2 i a . By 
using these distribution functions in Eq. (^) we can express the cross section as 

da da 



dQfdQ 1 dQ 2 h dQ 2 \dQ\dy x \dQ 2 dy 2 



(2vr) 8 1 1_ v 2 



X a lX a 2 XblXb 2 s al s a2 s bl s b2 



X 



J d 2 r 21a f a (x al , x a2 , r 2 i a , s al , s a2 )f b (x bl , x b2 , r 21a , s bl , s b2 ). (34) 
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If we assume that partons are weakly correlated and are evenly distributed in a hadron, we 
can use a flat distribution for the hadron wave function. For example, for hadron A we have: 



v a ,s a 



1 2 



Sa S a 



4x a ix a2 (2vr) 6 2 

-fix S -X S M(^-^) 2 -R 2 M(^-Ra) 2 -R 2 ) 
— Ja{ x al, Sal-, X a 2, S a2 ) {j[R 2 ) 2 ' 

By using this in Eq. ( p^) and by carrying out the d 2 v 2 i a integration, we get 
da da (2tt) 8 1 1 



dQtdQ^dQtdQ^ \dQ\dy 1 \dQ 2 dy 2 s 4 x al x a2 x bl x b2 SalS ^ blSh2 dp 7rR 2 

X f a {Xal,X a2 , S a i, S a2 )f b (x b i, X b2 , S bl , S b2 ), (36) 

where k w 0.5 is a geometrical factor associated with the overlapping area of the two colliding 
hadrons. For the cross section of the di-lepton production, we can use the result in Appendix 
A for \T\^ p and get 



da , ^ 2fi;e 2 e? „, . . w 87ra 2 87ra 



2 



- — 1 — — — \y 1=y2 =o = J2 ^^/Ki,z o2 )/(z&i,x 62 )(- )(- ), (37) 

dmidyidm 2 dy 2 ttK z 9m 1 s 9m 2 s 

where = \fQ] = sx a iX b i is the invariant mass square of the lepton pairs and e ai is the 
fractional charge of parton Oj. The summation is to sum over the flavors of parton a% and 
a 2 . The factor 2 in Eq. (|37j ) comes from when we include the Feynman diagram in Fig. 3 
and its cross diagram. Equation ( j37|) enables us to calculate the double parton cross section 
by using two universal double distribution functions and two individual " hard scattering" 
single parton cross sections 
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V. Gluon splitting process 



In this section we discuss the 4-final particle cross section from the gluon splitting 
process. We can use Eq. (|24] ) for this process if we take the interacting partons a\ and 61 as 
gluon partons. The corresponding partonic scattering amplitude square is shown in Fig. 4. 
There are two quark loop integrations in this Feynman diagram and one of them is 



d 2 \tdk + dk 



{2k+k~ - k 2 + ie)[2{k+ - Qt)(k~ - Q2 ) - (k - Q2) 2 + ie] 



1 

X [2(k+ + Qf)(k- - Qf) - (k - Q 2 ) 2 + ie]{2k+(k- - Q{ - Q 2 ) - k 2 + ie] ' (38) 

From the overall delta function, if we assume that all the transverse momenta of the partons 
are negligible, we have 

ai « (Qf + Qf, 0, 0), h « (0, + Q 2 , 0). (39) 



We first carry out the integration over k~ . When < k + < Q 2 , only the second 
k~-pole in the denominator is located on the lower half complex plane. We can close the 
contour on the lower plane and pick up a pole at 

* - « + W^-y (40) 

When —Qx < k + < 0, only the third A;~-pole in the denominator is located on the upper 
half complex plane. We can close the contour on the upper plane and pick up a pole at 

A plausible prejudgement is that the main contribution to the loop integration comes from 
when all the momenta in the loop are almost on their mass shell, i.e., when 

k- = Q 2 , k + = 0, k 2 = 0, (k-Q 2 ) 2 = 0. (42) 

With this in mind, we can set k 2 = k% = Q 2 and get 



27rzTr[ 7 ' 3 # A6 (^ - ^{qf + ft - (fc)£ Xa 



J-Qt (2A;+Q 2 --k 2 + 2e)2(Q+ + Q+)(k-Q 2 ) 2 (-2A;+Qr-k 2 + ze)' 1 1 

For Qf and Qf much larger than |k| and |k — Q 2 |, we can approximate the integration range 
of k + in Eq. ( f43"D to the whole space then integrate dk + by closing the contour in the upper 
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half plane and pick up a pole at k + = k 2 / (2Q 2 ). Eq. (H^) is then 
L = [ d 2 k 



±{Qt + Qt)(QT + )k 2 (k - Q 2 ) 2 

(44) 



There are many terms contribute to the denominator in Eq. fl44|) . One can show that the 
terms proportional to Q1Q1Q2Q2 are equal to zero. The leading contributions to the de- 
nominator in Eq. ([14]) are terms proportional to QfQ~k. l (k — Q2)"' and these terms make the 
integration over c? 2 k of Eq. ([0]) logarithmically divergent. One needs to introduce a cut-off 
to regulate this divergence. However, the cut-off is lost because of our previous approxima- 
tion. To mimic the physical cut-off we will use the (2 — 2e)-dimensional regularization in the 
MS -scheme. In dimensional regularization, one introduces a factor fi 2e to fix the dimension- 
ality. In order to simplify the fromula, one can express \i in terms of fi MS by the relation 
A*ms = 47r/i 2 e~ 7 B , in which 7^ is the Euler constant. The result of the d 2 ~ 2e h integration has 
a term with 1/e pole. The prescription is to drop this pole term, which effectively cuts off 
the intergration at |k| « fi MS . Since the physical cut-off is at k « Q 2 , we set /zms ~ | CJ 2 1 - 
Preceding with this prescription we can get 

d k F|k3Q^ = M( " ln ^ ) T-^ 1 ' (45) 



After some tedious calculations we get the leading contribution to the scattering matrix 
square integrated over the 4-lepton final momentum space as 

16a 2 a 4 7r 2 1 1 AQ 4 (iln 2 ^-ln^ + l) 
1 1 9(2vr)6 Q\Ql (Qt + Qt) 2 (Ql + Q2) 2 ' 

(46) 

We have inserted a color factor Cp = Tp/N g = 1/32 for Fig. 4 into Eq. (|46|), in which Tp = 
1/2 comes from 5 a pT F = Tr{t a tp) and iV s = 8 is the number of gluon color configurations. 
Inserting |T| 2 to Eq. (J2J) and assuming that for unpolarized hadrons the gluon distributions 
function of helicity 1 is equal to that of helicity -1, we get the cross section for the gluon 
splitting process as 

da da 
d 2 AQdQtdQidQtdQ2 = \&QdQ 2 d yi \dQ 2 2 dy 2 

= 2a 2 s a\2,f l 1 1 1 AQ 4 (lln 2 ^-ln^ + l) 

9s 2 x al x hl QlQ 2 2 {Ql + Qt) 2 {Qi+Q2) 2 

X fga(Xal)fgb(x bl ). (47) 
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VI. Discussions 



If we use the result in the Apendix for |T| 2 and integrate d 2 AQ for Eq. 



sp B ^. VC J) and Eq. 

(|i7|) up to a value AQ = q ma a; in which << q ma x << Qf, Qi , Qf > Q2 j we can see that 
the cross sections for single parton and gluon splitting processes behave like 

16e> 4 



da 

dm\dy\dm<idy2 



x 



E 



-fa(Xal)fb(Xbl)- 



7TQ 2 



3 9s 2 (27r) 2 x a ix w ' \jQ\yjQl 

(Qt + <&)(Qr + Q 2 ")(Qi-Qr + iQtQ~2 + Q^Q 2 
(QtQlQtQZ) 



(4$ 



da 



dniidy 1 dm2dy2 



^w /a(lal)/i(lM) ^^. 



MS; 



^(qL^/^Ms) + !] 



(^ + Q 2 f ) 2 (Qr+Q2) s 



(49) 



Since in these processes Qt ~ 0, AQ « Q x w — Q 2 , if we collect all the 4-particle 
events with transverse momenta Qi or Q 2 less than <\ maX) we should see the cross section of 
single parton process increases linearly with ^ max and the cross section of the gluon splitting 
process increases approximately as q^ m . However, from Eq. ( [37|) we see that the cross 
section of double parton process keeps at a constant value when f^ max increases. Also if we 
use the HERA result for the quark distribution function f(x) oc x~ 1A [§] and use Eq. (|l~5|) 
and Eq. (|26|) for the values of ar's, when we select the events with Qf , Qf, Qj, Qf ~Qwe can 
see that the differential cross section for single parton process hehaves like s 0A Q~ 6 ' 8 q L ^ nax and 
that for gluon splitting process behaves like s 0A Q~~ 1 ®' ? 'aft riax . For double parton cross section, 
if we assume that partons are evenly distributed in a hadron with radius R and assume 
that the double parton distribution function f(x\,X2) can be factorized into muliplication 
of two single parton distribution functions f(xi)f(x2), we see the differential cross section 
in Eq. (37) beheaves like s S Q~ 7s R~ 2 . For experiments of searching double parton events 
we suggest one do experiments at large colliding hadron center mass energy s and collect 
the events with small q^j. then check their cross section behavior on Q. For experiments of 
searching gluon splitting process we suggest one collect the events up to a large q max value 
and select the events with small Q. 

The significance of the result of the double parton cross sections in this paper are two- 
fold. We can use Eq. ([37]) for double parton cross section and fit it with the single parton 
distributions obtained from the deep inelastic scattering or Drell-Yan process then find the 
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hadronsize R or we can use the double parton experiments data to fit the new double parton 
distribution functions in a hadron. In both ways they give us a better understanding of the 
hadron structure. 
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APPENDIX 



In this appendix we calculate |T| 2 for single and double parton processes where |T| 2 is 
the partonic matrix square integrated over the 4-final particle phase spaces: 

4 r H 3 J- 

m 2 = (iy j^o)^(k + k-Qi)S% + h-Q 2 ).\m\ 2 (50) 
The following result is essential for calculating \T\ 2 : 

With this result, the leading Feynman diagram to |T| 2 p shown i Fig. 5 can be express as: 
1^,2 16e 4 a 4 QiaQia's, . Q2/3Q2/3' v 



X 



9(2tt)6 v * aa ' Q\ > y * pp ' Ql 
1 



Q?Q|(Qi-a) 2 (Qi-a) 2 
+ cross diagram contribution. (52) 

If we neglect the transverse momemta for partons, we can set a pa (Qf + $2,0, 0) and 
b pa (0, Qi + Q2 , 0). By using these in Eq. ([52] ) we get 



,2 lee 4 ^ 4 A(Qt + Qt)(Ql + Q 2 )(QjQi + 2QtQ 2 + QtQl 
lsp 9(2vr)6 QiQiQiQ 2 QlQl 



For double parton interaction, since it is a four-parton process in the feynman diagram, 
it is possible that the inteference terms from the spin flip processes as show in Fig. 5. may 
occur. The spin flip contribution as shown in Fig. 6 is: 

IT I 2 - — — Ql»Qla\ 

1 11 (2nfQj 3 1 9aa Ql ) 

x F*(Qr,Qi -a > -l/2)(-ie 9 7g 9 )u /9 (Q+ a, 1/2) 

x a'^/^^e^O^KQr^Qi-a',!/^. (54) 

To calculate 7\ we use spinors in Ref. for quarks moving in the +z direction, they are and 
use spinors in Eq. (|?p for quarks moving in the —z direction. The calculation showes that if 
we neglect the transverse momemta for partons, i.e , if we set a\ ~ (Qf, 0, 0), 61 ~ (0, Qi , 0), 
a 2 ~ (Q^,0, 0) and b 2 ~ (0,$2)0) the spin flip contribution is zero. With this result, we 
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can calculate |T| 2 for double parton process by using 

IT I 2 - ^4 1 271 1 , QiqQw x27T/ , Q2(iQ2(i' 

l dpl ~ (2vr) 12 g 2 gi 3 1 w+ g 2 ^ 3 l ^' + g 2 

x Tr[|5 l7 ^ l7 °']Tr^7^ 2 7' 3 '] 
+ cross diagram contribution. 



We get 



IT, 



dpi 



-2 "2 
, e al e a2 ( 



16 2 



(27T) 
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FIGURE CAPTIONS 



Fig. 1 Double parton interaction of hadron-hadron scattering. 

Fig. 2 Single parton interaction of hadron-hadron scattering. The dot line denotes quark 
or gluon partons. 

Fig. 3 Partonic scattering matrix square of double parton interaction. 

Fig. 4 Partonic scattering matrix square of gluon splitting process. 

Fig. 5 The leading Feynman diagram for single parton scattering. 

Fig. 6 Spin-flip process contribution to double parton interaction. 
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